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ABSTRACT 

The aim of the Degenerate Objects around Degenerate Objects (DODO) survey is to 
search for very low mass brown dwarfs and extrasolar planets in wide orbits around 
white dwarfs via direct imaging. The direct detection of such companions would allow 
the spectroscopic investigation of objects with temperatures much lower (< 500 K) 
than the coolest brown dwarfs currently observed. These ultra-low mass substellar 
objects would have spectral types >T8.5 and so could belong to the proposed Y dwarf 
spectral sequence. The detection of a planet around a white dwarf would prove that 
such objects can survive the final stages of stellar evolution and place constraints on 
the frequency of planetary systems around their progenitors (with masses between 
1.5 — 8 Mq, i.e., early B to mid F). This paper presents the results of a multi-epoch J 
band common proper motion survey of 23 nearby equatorial and northern hemisphere 
white dwarfs. We rule out the presence of any common proper motion companions, 
with limiting masses determined from the completeness limit of each observation, to 
18 white dwarfs. For the remaining five targets, the motion of the white dwarf is not 
sufficiently separated from the non-moving background objects in each field. These 
targets require additional observations to conclusively rule out the presence of any 
common proper motion companions. From our completeness limits, we tentatively 
suggest that < 5% of white dwarfs have substellar companions with T c g > 500 K 
between projected physical separations of 60 — 200 AU. 

Key words: stars: white dwarfs; planetary systems; low mass, brown dwarfs; imaging. 



1 INTRODUCTION 



Directly imaging the extrasolar planets found in orbit 
around solar type stars is difficult as these faint compan- 
ions are too close to their bright parent stars. As this pa- 
per was being finalised, Kalas et al. (2008) announced the 
discovery of a directly imaged ~ 3 Mj up extrasolar planet 
with a projected physical separation of 119 AU in orbit 
around the A- type star Fomalhaut. On the same day, Marois 
et al. (2008) announced the discovery of three directly im- 
aged companions around the A-type star HR8799 with likely 
masses between 5 — 13Mj up and projected physical sepa- 
rations of 24, 38 and 68 AU. However, coronagraphy and 
adaptive optics were needed to detect these faint extraso- 
lar planets. Another, perhaps simpler, solution to the prob- 
lems of contrast and resolution is to instead target intrin- 
sically faint stars. For example, many groups are already 



searching for planetary mass companions 1 in orbit around 
young, low mass stars and brown dwarfs (e.g., Chauvin et al. 
2003; Neuhauser et al. 2003). Any planetary mass com- 
panions found in orbit around these young stars will have 

1 We make the distinction between very low mass brown dwarfs 
and massive extrasolar planets by formation mechanism rather 
than mass, since the mass distributions of these two types of 
object likely overlap. For example, the 9 Mj up transiting object 
HAT-P-2b is too dense to be a brown dwarf (Baraffe et al. 2008), 
while the free floating objects with masses of the order of a few 
Jupiter masses that have been found in young clusters possibly 
formed in the same manner as stars. Indeed, some authors insist 
that the IAU distinction between these two populations, based on 
mass alone, has no valid foundation (Chabrier et al. 2008). There- 
fore, throughout this paper we prefer to use the term "planetary 
mass object" to refer to any body at or below the deuterium 
burning limit (13 — 14Afj up ), since the evolutionary history of 
any companions discovered with these masses is uncertain. 
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relatively high luminosities, since planetary mass objects 
cool continuously from the moment they form. Famously, 
a ~ 4 ± 1 Mjup (Ducourant et al. 2008) companion to the 
~ 25 Mjup brown dwarf member of the TW Hydrae associa- 
tion 2MASSW J1207334 - 393254 (2M1207) was imaged by 
Chauvin et al. (2004, 2005). However, Lodato et al. (2005) 
argue that 2M1207Ab more likely formed as a binary brown 
dwarf system, since the core accretion model, thought to be 
the most likely formation mechanism for gas giants like those 
in the Solar System, is unable to account for the formation 
of 2M1207b. 

An alternative approach, rather than looking at the 
bright, early part of a planet's life, is to look at the faint, 
late part of a star's life. White dwarfs are intrinsically faint 
stars and can be up to 10, 000 times less luminous than their 
main sequence progenitors, significantly enhancing the con- 
trast between any companion and the white dwarf. In ad- 
dition, any companion that avoids direct contact with the 
red giant envelope as the main sequence progenitor evolves 
into a white dwarf will migrate outwards as mass is lost 
from the central star by a maximum factor of Mms/Mwd 
(Jeans 1924). This increases the projected physical separa- 
tion between the companion and the white dwarf, substan- 
tially increasing the probability of obtaining a ground based 
direct image of a planetary mass companion. The evolution 
of planetary systems during the post-main sequence phase 
is discussed in more detail by Duncan & Lissauer (1998), 
Burleigh, Clarke & Hodgkin (2002), Debes & Sigurdsson 
(2002) and Villaver & Livio (2007). 

The direct detection of a planetary mass companion in 
orbit around a white dwarf would allow the spectroscopic in- 
vestigation of very low mass objects much cooler (< 500 K) 
than previously found. The coolest known brown dwarfs, 
ULAS J003402.77 - 005206.7 (Warren et al. 2007) and CF- 
BDS J005910.90 - 011401.3 (Delorme et al. 2008), have ef- 
fective temperatures of 600 < T e g < 700 K and spectral 
types of T8.5. The letter Y has been suggested for the next, 
cooler, spectral type (Kirkpatrick 2005) . Any planetary mass 
companions directly detectable around old (> 2 Gyr) white 
dwarfs could well belong to this class, regardless of forma- 
tion mechanism (Zuckerman & Song 2008). Such a discovery 
would help provide constraints on models for the evolution 
of planets and planetary systems during the final stages of 
stellar evolution. In addition, the age of any substellar and 
planetary mass companions discovered in such a system can 
be estimated using the white dwarf cooling age and the mass 
and the lifetime of the main sequence progenitor, providing 
model-free benchmark estimates of their mass and luminos- 
ity, which could be used to test evolutionary models (Pin- 
field et al. 2006). Radial velocity searches have concentrated 
mainly on stars with spectral types between mid F and M, 
since the faster rotation and increased activity of early B, A 
and mid F type stars broadens the low number of absorption 
lines in their spectra. As a result, it is difficult to accurately 
measure the Doppler shift of stars with these earlier spectral 
types. However, new methods in manipulating the measure- 
ments acquired when using the radial velocity technique has 
allowed planetary mass companions to be found around stars 
with spectral types of A and F (e.g., Galland et al. 2005). 
Nevertheless, as the 1.5 — 8Mq progenitor stars of white 
dwarfs have spectral types of early B, A and mid F, search- 
ing for planetary mass companions in orbit around white 



dwarfs allows the examination of a currently inadequately 
explored region of parameter space, supplying new informa- 
tion on the frequency and mass distribution of extrasolar 
planets around intermediate mass main sequence stars. 

A number of extrasolar planets have been discovered 
around evolved giant stars using the radial velocity tech- 
nique, e.g., HD 11977 (G5 III; Setiawan et al. 2005), 
HD 13189 (K2 II; Hatzes et al. 2005) and /3 Gem (K0 III; 
Hatzes et al. 2006; Reffert et al. 2006). These stars have 
entered the red giant phase of stellar evolution, proving 
that planets can survive the early stages of the RGB phase. 
Evolved giant stars are significantly more massive than the 
Sun, so their progenitors were likely to be A or B type stars 
(see Table 6 of Hatzes et al. 2006) . The companions all have 
masses significantly larger than Jupiter, implying that sig- 
nificantly more massive planets are formed around these in- 
termediate mass stars than around lower mass stars (Lovis 
& Mayor 2007). In fact, both Lovis & Mayor (2007) and 
Johnson et al. (2007) suggest that intermediate mass stars 
are more likely to host extrasolar planets of all masses than 
solar mass stars. 

Up to three white dwarfs are known to be wide com- 
panions to stars hosting extrasolar planets. The bright 
(V = 11 mag), well studied white dwarf WD 1620 - 391 
(CD— 38° 10980) was known to be a common proper motion 
companion to the solar type star HD 147513 (Wegner 1973) 
before a planet, with a minimum mass, M p sin i = 1.21 Mj up 
and an orbital radius of 1.32 AU, was discovered in orbit 
around the latter (Mayor et al. 2004). Since WD 1620-391 
and the parent star are separated by ~ 5360 AU, it is highly 
unlikely that the evolution of the main sequence progen- 
itor of the white dwarf affected the mass or the orbit of 
HD 147513b. It has been suggested that the faint com- 
panions to two other planet hosting stars are also white 
dwarfs. The planet in orbit around Gliese 86 has a mini- 
mum mass, MpSini = 4.01 Mj up and an orbital period of 
15.8 days (Queloz et al. 2000), while the likely white dwarf 
companion has a mass between 0.48 5j M ^ 0.62 Mq and 
an orbital radius of ~ 20 AU (Mugrauer & Neuhauser 2005; 
Lagrange et al. 2006). The possible effects of the evolution 
of the main sequence progenitor on the mass and the orbit 
of the planet are discussed by Desidera & Barbieri (2007). 
The planet in orbit around HD 27442 (e Ret) has a min- 
imum mass, MpSini = 1.28 Mj up and an orbital radius of 
1.18 AU (Butler et al. 2001). The white dwarf companion 
(Raghavan et al. 2006; Chauvin et al. 2006), currently sepa- 
rated from HD 27442 by ~ 240 AU, was recently confirmed 
from an analysis of its optical and infrared (IR) spectrum 
(Mugrauer et al. 2007; Chauvin et al. 2007). 

The discovery of metal rich dust disks in close orbits 
around white dwarfs may indicate the existence of old, rocky 
planetary systems, suggesting that even terrestrial planets 
and asteroids can survive the final stages of stellar evolu- 
tion. The first dust disk was discovered around the DAZ 
white dwarf G 29 - 38 (WD 2326 + 049) from the identifi- 
cation of a large IR excess in its spectrum at wavelengths 
between 2 — 5/im (Zuckerman & Becklin 1987b). This IR 
excess was initially attributed to a spatially unresolved, 
T e ff = 1200 ± 200 K brown dwarf companion to the white 
dwarf. However, subsequent measurements showed that the 
blackbody-like IR excess was due to a dust disk rather than 
a brown dwarf (Tokunaga, Becklin & Zuckerman 1990). A 
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mid-IR (MIR) spectrum of G 29 - 38, obtained by the 
Spitzer Space Telescope, shows a strong emission feature in 
the spectrum between 9 — 11/im, which indicates the pres- 
ence of silicates (SiO.4) in the dust disk (Reach et al. 2005). 
Eight additional white dwarfs are now known to have dust 
disks in orbit around them: GD 362 (T cff = 9740 K; Beck- 
lin et al. 2005; Kilic et al. 2005), GD 56 (T cff = 14400 K; 
Kilic et al. 2006), WD 1150 - 153 (T cS = 12800 K; Kilic 
& Redfield 2007), WD 2115 - 560 (T off = 9700 K; von 
Hippel et al. 2007), GD 40 (T eff = 15200 K; Jura, Farihi 
& Zuckerman 2007), GD 133 (T eff = 12200 K; Jura et al. 
2007), PG 1015+161 (Tea = 19300 K; Jura et al. 2007) and 
G 166-58 (T eff = 7390 K; Farihi et al. 2008). The generally 
accepted origin of the material in these dust disks is from 
the tidal disruption of an asteroid that had strayed within 
the Roche lobe radius of the white dwarf after its orbital ra- 
dius was altered during the AGB phase of stellar evolution 
(Graham et al. 1990; Debes & Sigurdsson 2002; Jura 2003). 
In addition to these dust disks, there have been metal rich 
gas disks found in orbit around 2 hotter DAZ white dwarfs; 
SDSS J122859.93 + 104032.0 (T eff = 22292 K; Gansicke et al. 
2006) and SDSS J104341.53 + 085558.2 (T eff = 18330 K; 
Gansicke, Marsh & Southworth 2007). These gas disks could 
also indicate the presence of old planetary systems, since it 
is likely that the hot white dwarfs caused the dust in the disk 
to sublimate. The fraction of known single DAZ white dwarfs 
with IR excesses, which can be attributed to a dust disk, is 
14% (Kilic et al. 2006). In addition, Jura (2006) argues that 
> 7% of white dwarfs possess asteroid belts similar to that 
of the Solar System, and by implication, remnant planetary 
systems. 

The first search for low mass substellar companions to 
white dwarfs was conducted by Probst (1983), who mea- 
sured the IR magnitudes of ~ 100 white dwarfs to determine 
whether any excess emission was present. No companions 
were found during this survey. Subsequently, a number of 
groups unsuccessfully attempted to detect substellar com- 
panions to white dwarfs using the same method (e.g., Ship- 
man 1986; Zuckerman & Becklin 1987a). The first confirmed 
subtellar companion to a white dwarf was discovered in 1988 
around the DA white dwarf GD 165 (Becklin & Zuckerman 
1988). Over 15 years later, a second brown dwarf compan- 
ion was found in orbit around GD 1400 (Farihi & Christo- 
pher 2004; Dobbie et al. 2005). More recently, radial veloc- 
ity measurements revealed a brown dwarf companion in a 
close (~ 116 minutes) orbit around WD 0137-349 (Maxted 
et al. 2006), while its spectral type (L8) was later determined 
from near-IR (NIR) spectroscopy (Burleigh et al. 2006). The 
L dwarf companion fraction, determined from a wide field, 
proper motion, NIR search for wide substellar companions 
to 347 white dwarfs, is estimated to be < 0.5% (Farihi et al. 
2005). 

The recent discovery of an extrasolar planet around 
the post-red giant star V 391 Pegasi proves that plane- 
tary mass companions with an initial orbital radius outside 
the maximum radius of the red giant envelope can survive 
the RGB phase of stellar evolution. The planet was discov- 
ered from the periodic variation in the precise timing mea- 
surements of V 391 Pegasi's extremely stable, short period 
pulsations (Silvotti et al. 2007). It has a minimum mass, 
MpSini ~ 3.2Mj up , an orbital radius of ~ 1.7 AU and an 
estimated age of ~ 10 Gyr. Strong evidence for the existence 



of a planetary mass companion to the DAV white dwarf 
GD 66 has been recently found from the periodic variation 
in the precise timing measurements of GD 66's extremely 
stable non-radial pulsations (Mullally et al. 2008). While 
current measurements suggest that this companion has a 
minimum mass, M p sini ~ 2.11 Mj up and an orbital radius 
of ~ 2.356 AU, further observations, to cover the entire orbit 
of the companion, are now required. 

In Burleigh et al. (2008) we reported limits on planetary 
mass companions to the nearest single white dwarf, vMa 2. 
Preliminary results and progress reports from the Degen- 
erate Objects around Degenerate Objects (DODO) survey 
have been published previously (Clarke & Burleigh 2004; 
Burleigh, Hogan & Clarke 2006; Hogan, Burleigh & Clarke 
2007). In this paper we report further extensive results of 
the DODO survey; a NIR direct imaging search for substel- 
lar and planetary mass common proper motion companions 
in wide orbits around nearby white dwarfs. 



2 TARGET SELECTION 

The ability to directly image an extrasolar planet in orbit 
around a white dwarf will depend on the apparent magni- 
tude of the planet, which in turn depends on its absolute 
magnitude and distance from the Earth. The absolute mag- 
nitude of the planet is determined from its intrinsic luminos- 
ity, which is dependant primarily on the age and the mass of 
the planet, since it will cool continuously from the moment 
it formed. The age of an extrasolar planet found in orbit 
around a white dwarf equals the sum of the main sequence 
progenitor lifetime and the white dwarf cooling age. Using 
this age and the distance to the white dwarf, Burleigh et al. 
(2002) used evolutionary models for cool brown dwarfs and 
extrasolar planets (Burrows et al. 1997) to make initial pre- 
dictions of the IR magnitudes of planetary mass companions 
around white dwarfs. The results published in this paper use 
the more recent "COND" evolutionary models of Baraffc 
et al. (2003). These models assume irradiation effects from 
the parent star on the planet are negligible and predict that 
a 5 Mjup planet with an age of ~ 2 Gyr will have an appar- 
ent magnitude of J ~ 24 mag at 10 pc. This magnitude is 
comparable with the expected sensitivity of a one hour ex- 
posure acquired using an 8m telescope. Brighter companions 
such as massive brown dwarfs and M dwarfs should easily be 
detected. Indeed, these more luminous objects would have 
already been detected in previous IR studies (e.g., Farihi 
et al. 2005) and in 2MASS data. 

Lower mass planets could be more easily detected 
around younger white dwarfs, since these companions would 
be brighter than their older counterparts. In addition, 
fainter, and therefore lower mass, companions could be more 
easily detected around nearby white dwarfs compared to 
more distant targets. Nearby white dwarfs are also more 
likely to have large proper motions, which require a smaller 
baseline between the observations of the two epochs. An ini- 
tial sample of ~ 40 targets, with total ages (main sequence 
progenitor lifetime plus the white dwarf cooling age) < 4 Gyr 
and distances within ~ 20 pc, were selected from the cata- 
logue of white dwarfs within 20pc Holberg et al. (2002). Of 
these targets, 23 equatorial and northern hemisphere white 
dwarfs (Table 1) are presented in this paper. The remaining 
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Columns: fj, and 9 are the R.A. and Dec components of the proper motion of the white dwarf, respectively, measured in milli arc 
seconds per year; d is the distance to the white dwarf, measured in parsecs; T e ff is the effective temperature of the white dwarf, 
measured in Kelvin; log g is the log of the gravity of the white dwarf; Mwd is the mass of the white dwarf, measured in solar 
masses; %d is the cooling age of the white dwarf, measured in gigaycars; Mms is the mass of the main sequence progenitor, 
measured in solar masses, and is calculated using the IFMR of Dobbic ct al. 2006; £ms is the main sequence lifetime, measured in 
gigaycars, and is calculated using Equation 2 (Wood 1992); ttot is the total age of the white dwarf, measured in gigayears. 
R = References, which refer to the T ef j, log g, M^q and columns. (1) Salim & Gould (2003), (2) van Altcna, Lee & HofHcit 
(1995), (3) Dufour, Bergeron & Fontaine (2005), (4) Perryman et al. (1997), (5) Bergeron, Saffer & Liebert (1992), (6) Bergeron, 
Leggett & Ruiz (2001), (7) Fontaine, Bergeron & Brassard (2007), (8) Liebert, Bergeron & Holbcrg (2005), (9) Bakos, Sahu & 
Nemeth (2002), (10) Gianninas, Bergeron & Fontaine (2005), (11) Holbcrg, Oswalt & Sion (2002), (12) Koester et al. (2001), (13) 
Giovannini ct al. (1998), "Magnetic field strength, B = 100 MG; rotation period, P = 1.3 days (Wickramasinghe & Ferrario 2000), 
^Magnetic field strength, B = 320 MG (Wickramasinghe & Ferrario 2000), c Magnetic field strength, B = 70 kG; rotation period, 
P = 1.4418 days. WD 1953 — 011 is also photometrically variable at the ~ 2% level, an effect which is believed to be caused by 
a star spot (Brinkworth et al. 2005), d A dust disk is known to orbit this ZZ ceti white dwarf, *The white dwarf cooling age was 
calculated using models from Fontaine, Brassard & Bergeron (2001). 

southern hemisphere white dwarfs will be presented in an 
upcoming paper. 

The cooling age of a white dwarf, iwD, can be calcu- 
lated using evolutionary models. When the cooling age was 
unavailable in the literature, models from Fontaine, Bras- 
sard & Bergeron (2001), which use T e ff and log g values to 
calculate the cooling age, were used to estimate this value. 
The initial final mass relation (IFMR) determined by Dobbie 
et al. (2006), based on the measurements of a small number 
of white dwarfs found in young open clusters, was used to 
determine the mass of the main sequence progenitor, Mms, 
from the mass of the white dwarf, Mwd- This linear IFMR 
is given as 

Mwd = 0.133 Mms +0.289 (1) 

Recent observations of white dwarfs in older open clusters 
have placed constraints on the low mass end of the IFMR, 
suggesting that this equation is valid down to white dwarf 
masses of 0.54 M (Kalirai et al. 2008). The main sequence 
progenitor lifetime, £ms , is estimated using the equation 



where tyis is measured in Gyr (Wood 1992). 



3 OBSERVATIONS 

Observations of the 23 equatorial and northern hemisphere 
white dwarfs (Table 2) were acquired in the J band primarily 
using Gemini North and NIRI between 2003 and 2005, while 
a small number of observations of equatorial targets were 
acquired in 2002, using Gemini South and FLAMINGOS. 
NIRI consists of a 1024 x 1024 pixel ALADDIN-II array. 
When combined with the f/6 camera, NIRI supplies a pixel 
scale of 0.117" pixel" 1 and a wide field of view of 120" x 120". 
FLAMINGOS consists of a 2048 x 2048 pixel HAWAII-II 
array. When combined with the f/16 camera, FLAMINGOS 
supplies a pixel scale of 0.078" pixel -1 and a wide field of 
view of 160" x 160". 
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Table 2. Details of the observations of the 23 equatorial and 
northern hemisphere white dwarfs in the DODO survey 



White Dwarf Date ET FWHM TW 

Number Observed [m] ["] 



The science data were acquired using a dither pattern, 
which involved systematically offsetting the telescope, to al- 
low the effective removal of the sky background. The total 
exposure time given in Table 2 was achieved by obtaining 
60 second and 90 second individual exposures per dither po- 
sition for NIRI and FLAMINGOS, respectively. The number 
of coadds acquired for each individual exposure was adjusted 
to avoid saturating the white dwarf. Unfortunately, on occa- 
sion, the requested exposure time of 1 hour was not always 
achieved. "Lamps on" dome flats were acquired by imaging 
a uniformly illuminated screen within the dome. "Lamps 
off" dome fiats were acquired using the same method, ex- 
cept with no illumination of the screen. These dome flats 
were used for the calibration of the NIRI science data. High 
and low twilight flats were acquired for the calibration of the 
FLAMINGOS data. Short dark frames were also acquired 
to help with the identification of bad pixels. 



4 DATA REDUCTION 

All the data acquired for the DODO survey were reduced 
using the Image Reduction and Analysis Facility (IRAF; 
Tody 1986) and the Gemini package, versions 2.12.2a and 
1.8, respectively. Raw NIRI images are in the form of multi- 
extension fits (MEF) files with most of the header informa- 
tion in the Primary Header Unit (PHU), "[0]" extension 
and the raw image data in the second, "[1]" extension. Raw 
FLAMINGOS images are in the form of single extension 
fits files. The nprepare and fprepare tasks in the Gemini 
package were applied to all raw data acquired with NIRI 
and FLAMINGOS, respectively. These tasks add certain es- 
sential keywords to the header of each data file, allowing the 
subsequent data reduction tasks to be applied. In addition, 
the fprepare task converted all the FLAMINGOS images 
into MEF files. This made it possible for both the NIRI and 
the FLAMINGOS data to be reduced in a homogeneous 
manner, using the niri tasks in the Gemini package. 

For the reduction of the science data, a sky frame was 
used instead of a dark frame as it gives a much better dark 
measurement, since it was acquired concurrently with the 
science data. The sky frame also removes any constant ad- 
ditions to the science data due to bias levels. The nisky 
task was used to create the sky frame by median combin- 
ing each individual science image after masking the objects 
in each image. The nireduce task was used to subtract the 
sky frame from the individual science images. A flat field im- 
age was created by subtracting the median combined "lamps 
off" / low twilight flat image from each individual "lamps 
on" / high twilight flat images and then median combing 
the resultant images. This flat field image was divided by 
its mean pixel value to create a normalised flat field im- 
age. The niflat and nireduce tasks were used to create 
this normalised flat field image and to divide the individual 
science images by the normalised flat field image, respec- 
tively. For each individual science image, a sky background 
image was created by median combining the previous and 
subsequent 5 and 4 science images for the NIRI and the 
FLAMINGOS data, respectively, after masking the objects 
in each image. The XMOSAIC task in the xdimsum package 
was used to create and subtract each sky background im- 
age from each individual science image and then create an 
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Columns: ET is the total exposure time, measured in min- 
utes; FWHM is the full-width at half-maximum, calculated 
as the average FWHM of stars in the field using SExtrac- 
tor, measured in arc seconds; TW is the telescope and in- 
strument the data were taken with; GN+N indicates Gemini 
North and NIRI were used; GS+F indicates Gemini South 
and FLAMINGOS were used; Notes: (1) Moderate 60Hz sig- 
nal, (2) Severe 60Hz signal, (3) Streak across the image due 
to a bright source just outside the field of view. 
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average combined final stacked image. This task also cor- 
rected the bad pixels in each sky subtracted image by using 
a mask created by the niflat task and linearly interpolating 
across bad columns in each image line. In addition, cosmic 
ray events are removed by replacing the value of the pixels, 
which are more than 5<r above the median of the pixel values 
in a surrounding box of 5 x 5 pixels and that are not part 
of an object, with the local median value. 

Saturated pixels (e.g., from the cores of bright stars) can 
lead to persistence effects in subsequent frames as residual 
charge remains in the pixels after the array has been reset. 
Persistence manifests itself as an apparent faint object at 
the location of the saturated object in the previous image. 
This faint object can remain for several frames depending 
on the level of saturation. For the ALADDIN-II NIRI array, 
a highly saturated star will leave a faint object at the level 
of ~ 1% in the subsequent exposure and ~ 0.2% in the next 
one (Hodapp et al. 2003). Persistence affects ~ 42% of the 
science data in Table 2. To remove these persistence effects 
from the sky background image and from the final stacked 
image, a mask consisting of the cores of the bright stars in 
the field from the previous three science images was created 
for each individual science image and added to the object 
mask used in the XMOSAIC task. 

Intermittent pattern noise substantially degraded a 
large amount of the NIRI data throughout 2003 and has 
also affected some of the data from 2004 and 2005. A diago- 
nal herringbone pattern due to 60 Hz interference (Hodapp 
et al. 2003) is seen in ~ 37% of the NIRI data (Table 2). This 
pattern is reflected symmetrically in each quadrant, due to 
the fact that the readout of the quadrants, from each corner 
of the array to the centre of the array along each row, is 
symmetrical. This pattern is not present in the lab and has 
been eliminated at times on the telescope, indicating that 
it arises from the telescope environment and not from the 
NIRI electronics. It was not possible to remove this noise 
from the NIRI data. 

Another form of intermittent pattern noise found in the 
NIRI data is due to 50 Hz interference, which creates a 
horizontal pattern. The constant variations in the rows and 
columns of each image due to this pattern were removed 
by collapsing each quadrant along rows and columns and 
then subtracting the median value from each row and each 
column separately for each of the quadrants. 

Another form of intermittent pattern noise found in the 
NIRI data, which manifests itself as a combined quadrant 
pattern plus a vertical pattern, was also removed in the same 
manner. The quadrant pattern occurs due to mismatched 
bias levels between the quadrants and the vertical pattern 
has a periodicity of 8 pixels, which corresponds to the 8 
amplifier channels reading out each quadrant. 



5 DATA ANALYSIS 

The cleaned final stacked images created using XMOSAIC 
were astrometrically calibrated to within ~ 1 — 2" using 
Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006) 
objects in the field of each white dwarf. The IRAF tasks 
CCFIND and CCMAP in the IMAGES package were used to cal- 
culate the transformation required to match the positions of 
the objects in the final stacked image with the positions of 



the objects in the 2MASS catalogue. Objects near the edge 
of the final stacked image were removed from the transfor- 
mation calculations. The ccsetwcs task was then used to 
apply this transformation to the final stacked image. The as- 
trometrically calibrated final stacked images were then pho- 
tometrically calibrated using aperture photometry. Objects 
in the final stacked image, with instrumental magnitudes 
mi, were matched with objects in the 2MASS catalogue. 
The instrumental magnitude of the matched object could 
then be associated with their apparent magnitudes, m, us- 
ing the linear formula m = m\ + zp. The zeropoint, zp, is 
therefore equal to the y intercept of the line of best fit to 
the points in the plot of m\ against m. Only photometric 
2MASS stars, with a J band photometric quality flag equal 
to "A" (given to objects with a SNR > 10 and a photomet- 
ric uncertainty, a < 0.109), were used to determine the zp. 
In addition, objects near the edge of the final stacked image 
or that were saturated were excluded from this step. 



5.1 Point Source Detection 

The SExtractor program (Bertin & Arnouts 1996) was used 
to detect all objects in the final stacked images with a signal 
to noise ratio (SNR) > 3. SExtractor determines the posi- 
tion of the centre of an object by computing the first order 
moments of the isophotal profile of the object, which is ad- 
equate for detecting point sources in all the final stacked 
images in the DODO survey. A "weight map", created by 
the XMOSAIC task, was used to normalise the noise over the 
final stacked image. This avoided the detection of spurious 
sources around the edges of the final stacked image, where 
only a few of the individual science images contribute. The 
first pass of SExtractor used apertures with diameters rang- 
ing from 1 to 20 pixels to determine the aperture size that 
delivered the highest SNR. Only objects with an internal 
flag equal to 0, indicating that no problems were found with 
the detection, were used in the determination of the opti- 
mum aperture size. In addition, the ellipticity of the objects, 
given as 1 — B/A, where A and B are the semi-major and 
semi-minor axes of the object, respectively, was chosen to 
be < 0.2. This excluded objects with high ellipticities, such 
as background galaxies, from the determination of the opti- 
mum aperture size. Also, only those objects within the cen- 
tral region of the final stacked image that have a full-width 
at half-maximum (FWHM) < 1" were used. This further 
assisted in removing extended objects. The resulting ideal 
aperture size was then used to detect all objects in the final 
stacked images. 



5.2 Measurement of Proper Motions 



The motion of the objects in the field of each white dwarf 
between the first epoch and second (or third) epoch images 
was calculated. Since the white dwarf is rarely positioned 
on the same pixels in each epoch, spurious distortion effects 
can be seen, which are caused by optical aberrations. As a 
result, the motion of objects between the two epoch images 
is seen to be a function of field position (Figure 1). When 
a large number of background reference stars are present, 
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Figure 1. Distortion effects before (left) and after (right) the distortion correction was applied to the images of WD 2047 + 372. The 
x- and y-axes are the R.A. and Dec, respectively, measured in degrees. The arrows in the image indicate the direction and magnitude, 
multiplied by 20, of the motion of each object in the field between the first epoch and second epoch images. 



the two epoch images can be well matched and the distor- 
tions can be effectively removed. However, as the number 
of background reference stars decreases, the two epoch im- 
ages can not be accurately matched and this is often the 
limiting factor for astrometric accuracy. The GEOMAP and 
GEOXYTRAN tasks in the IMAGES package were used to cor- 
rect for these distortion effects. Objects from the first epoch 
image were matched to the closest object present in the sec- 
ond (or third) epoch image only if 1) their magnitudes are 
within 1 mag, 2) the SExtractor internal flag ^ 3, which 
indicates a good detection, 3) the ellipticity of the object 
is < 0.5. This excluded objects with very high ellipticities, 
such as background galaxies, from the matching procedure. 
In some cases, the closest object was too far away to be a 
true match, so a clipping factor was introduced to remove 
these mismatches. 



5.3 Limits and Errors 

The completeness limit for each final stacked image was es- 
timated by determining the magnitude at which 90% and 
50% of inserted artificial stars were recovered from each 
image. The STARLIST task was used to create a list of 
200 randomly positioned artificial stars at a magnitude of 
J = 19.0 mag. The mkobjects task was used to insert 
the artificial stars into the final stacked image. SExtractor 
was then used to detect all objects in the image, includ- 
ing the artificial stars. The calculated magnitudes of the 
artificial stars were checked to ensure they were equal to 
J = 19.0 mag. Using the same artificial star list, the mkob- 
jects and SExtractor steps were repeated for magnitudes 
between 19.1 ^ J ^ 24.0 mag in 0.1 magnitude steps. The 
entire process was then repeated a further 50 times, equiva- 
lent to a total of 10, 000 inserted artificial stars for each 0.1 
magnitude bin. Plots of the percentage of artificial stars re- 
covered against the apparent J magnitudes of the artificial 
stars were created. The number of artificial stars recovered 
was often much less than 100% at the brighter J magnitudes 
as some stars were lost within the point spread function 



(PSF) of other real objects or artificial stars. The motion of 
an object can be calculated only when the object is detected 
in both epochs. Assuming that the probability of detecting 
an object in the first epoch image, Pi, is independent from 
the probability of detecting an object in the second epoch 
image, P2, the probability of detecting an object in both 
epochs is Pi x P2 . Therefore, by multiplying the individual 
completeness limits for each epoch, a combined complete- 
ness limit for both epoch images can be determined. This 
assumption is valid for objects near the completeness limit. 
However, it is not valid for the bright objects not detected 
due to the fact that they are within the PSF of other real 
objects or artificial stars. Therefore, the completeness limits 
at the brighter J magnitudes are underestimated. 

The "COND" evolutionary models for cool brown 
dwarfs and extrasolar planets (Baraffe et al. 2003), along 
with the magnitudes at which 90% and 50% of artificial 
stars were recovered, were used to estimate the minimum 
mass of a companion which could be detected in both epoch 
images. The models predict the absolute magnitudes of sub- 
stellar objects depending upon their age. Isabelle Baraffe 
kindly supplied these models for the total ages determined 
for all the white dwarfs in the DODO survey. The total age 
is equal to the sum of the main sequence progenitor lifetime 
and white dwarf cooling age, both of which depend upon 
evolutionary models. While the cooling age errors are small 
and well constrained (Fontaine et al. 2001), and the scatter 
in the empirical IFMR is significantly reducing as more and 
higher quality observations are made of white dwarfs in open 
clusters (Casewell et al. 2008) , the main sequence progenitor 
lifetimes rely on models which are difficult to calibrate (e.g., 
Catalan et al. 2008). Therefore, to take these uncertainties 
into account, a conservative error of ±25% is applied to the 
total age of each white dwarf (note that the white dwarf 
cooling age is the dominant timescale for most of the tar- 
gets in the DODO survey, as shown in Table 1). However, 
at ages > 1 Gyr, the "COND" evolutionary models indi- 
cate that the absolute magnitudes of substellar objects are 
relatively insensitive to changes in their age, implying that 
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even with a ±25% error, the resulting error on the mass of 
a companion is small (Table 3). 

The detection of a companion with a mass equal to 
the minimum mass determined using the completeness limit 
will only be possible if the companion is outside the extent 
of the PSF of the white dwarf. In addition, it is expected 
that the orbital radius of any companions that avoid di- 
rect contact with the red giant envelope will expand, which 
would increase the projected physical separation between 
the companion and the white dwarf. The majority of the 
DODO survey observations were acquired in good seeing 
conditions, so the minimum projected physical separation at 
which a companion could be found around each white dwarf 
was taken to be 3". Beyond this distance, the contribution of 
the flux from the white dwarf was assumed to be minimal. A 
more careful treatment of the PSF of the white dwarf could 
allow companions to be uncovered within 3" and will be 
dealt with in a future publication. The maximum projected 
physical separation at which a companion could be found 
around each white dwarf is limited by the field of view cov- 
ered by both epochs. The completeness limit is only valid 
in the central region of each final stacked image, where all 
the individual images contribute. The useable field of view 
decreases further when the two epoch images are matched 
as the white dwarf is rarely positioned on the same pixels in 
each epoch image. The minimum and maximum projected 
physical separations at which a companion could be found 
around the main sequence progenitor can also be estimated, 
since the orbital radius of any companions around the main 
sequence progenitor will expand by a factor of Mms/Mwd 
during stellar evolution. 



6 RESULTS 

Figures 2-24 show the results for each of the 23 equatorial 
and northern hemisphere white dwarfs from the DODO sur- 
vey. Using the combined completeness limit for each white 
dwarf, an estimate of the minimum mass of a companion 
which could be detected in both epoch images is calculated. 
The range of projected physical separations at which a com- 
panion of this mass could be found around each white dwarf 
and the corresponding range of projected physical separa- 
tions around the main sequence progenitors is determined. 
These results are summarised in Table 3. Comments arc 
made only on interesting objects. 

6.1 WD 0644 + 375 

The mass of WD 0644 + 375 used throughout this paper was 
determined by assuming that the core of this white dwarf is 
made partly of strange matter (Mathews et al. 2006). This 
unusual core composition was suggested as a way to explain 
the inconsistency between the radius determined from the 
parallax of WD 0644 + 375, obtained from Hipparcos data, 
and the radius predicted using a mass-radius relation, which 
assumes the core of the white dwarf is composed primarily 
of carbon (Provencal et al. 1998). The mass of the white 
dwarf was originally determined by Bergeron et al. (1992) 
to be 0.66 Mq, but this predicted a radius that was signif- 
icantly larger than predicted using the parallax. Therefore, 
the slightly lower mass of 0.54 Mq determined by Fontaine 



et al. (2007) is used here to determine the total age of the 
white dwarf, since this mass provides a radius consistent 
with observations. Note that a companion with a mass of 
5±1 Mj U p could have been detected if the larger white dwarf 
mass was used to determine the total age of the white dwarf. 

6.2 WD 0738 - 172 

WD 0738 — 172 is a member of a known common proper 
motion binary. The secondary star of this binary system is 
an M6 main sequence star (Monteiro et al. 2006) with an 
orbital radius of ~ 262 AU (Poveda et al. 1994). The main 
sequence secondary does not appear in the proper motion 
diagram as it was saturated in the 2005 second epoch image, 
making it unavailable for proper motion measurements. The 
overall decrease in the completeness limit, compared to the 
other white dwarfs, of the images acquired of WD 0738 — 172 
is due to the higher proportion of artificial stars inserted 
within the PSF of the bright secondary. 

6.3 WD 1626 + 368 

Recent MIR observations of the helium atmosphere DZ 
white dwarf WD 1626 + 368 show no evidence of a dust disk 
(Mullally et al. 2007). However, the abundance of carbon 
relative to iron in the atmosphere of WD 1626 + 368 is 10 
times below the solar abundance, similar to the carbon defi- 
cient asteroids in the Solar System. Therefore, external pol- 
lution from such asteroids naturally explains the abundances 
of the metals in the atmosphere of this white dwarf (Jura 
2006). The possible presence of asteroids in orbit around 
WD 1626 + 368 represents an increased probability of the 
existence of an old planetary system. However, the motion of 
this white dwarf between the first epoch and second epoch 
images is large enough to confidently state that there are 
no common proper motion companions to WD 1626 + 368 
within the limits given in Table 3. 

6.4 WD 1633 + 433 

Although no dust disk has been found in orbit around the 
DAZ white dwarf WD 1633 + 433, the presence of metals in 
its atmosphere may indicate the existence of an old plane- 
tary system. The 2003 first epoch image of WD 1633 + 433 
was degraded by 60 Hz interference (Table 2) , which has de- 
creased the completeness limit of this image. In addition, the 
2004 second epoch image of WD 1633 + 433 was degraded 
by a large streak across the image, due to a source just out- 
side the field of view. However, this streak is not present 
in the first epoch image. It is likely that these effects have 
introduced the large scatter in the motions of objects with 
magnitudes J > 21 mag between the 2003 first epoch and 
2004 second epoch images (Figure 25), particularly in the 
region of the streak. This suggests that the error on the 
motion of these faint objects is comparable to the motion 
of WD 1633 + 433 (Table 1). As a result, multiple objects 
appear to have motions similar to the motion of the white 
dwarf (the two objects with motions closest to the motion of 
WD 1633 + 433 lie on the streak). Real common proper mo- 
tion companions to WD 1633 + 433 cannot be distinguished 
from non-moving background objects. Therefore, a third 
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Completeness limits for the images of WD0115+159 Proper Motion Diagram for WD0115+159 
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Figure 2. The completeness limit (left) shows the percentage of artificial stars recovered by SExtractor from the 2003 first epoch and 
2004 second epoch images acquired of WD 0115 + 159, against the apparent J magnitude of the artificial stars. The proper motion 
diagram (right) shows the motion of all objects in the field of WD 0115 + 159 between the first epoch and second epoch images. The 
dashed green circles represent the la and 3cr scatter of the distribution of the motions of all objects excluding the white dwarf, centred 
on the white dwarf, to help determine possible common proper motion companions to the white dwarf. 



Completeness limits for the images of WD0148+467 Proper Motion Diagram for WD0148+467 
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Figure 3. The completeness limit (left) and the proper motion diagram (right) for WD 0148 + 467. 
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Completeness limits for the images of WD0208+396 Proper Motion Diagram for WD0208 + 396 
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Figure 4. The completeness limit (left) and the proper motion diagram (right) for WD 0208 + 396. 
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Completeness limits for the images of WD0738— 172 
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Figure 8. The completeness limit (left) and the proper motion diagram (right) for WD 0738 — 172. 




Figure 9. The completeness limit (left) and the proper motion diagram (right) for WD 0912 + 536. 
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Completeness limits for the images of WD1344+106 Proper Motion Diagram for WD1344+106 



100 - 




19 20 21 22 23 24 -500 

Apparent J Magnitude (mag) 5RA (mas) 

Figure 13. The completeness limit (left) and the proper motion diagram (right) for WD 1344 + 106. 



Completeness limits for the images of WD1609+135 Proper Motion Diagram for WD1609+135 



100 - 




19 20 21 22 23 24 -200 200 4C 

Apparent J Magnitude (mag) 5RA (mas) 

Figure 14. The completeness limit (left) and the proper motion diagram (right) for WD 1609 + 135. 



Completeness limits for the images of WD1626+368 Proper Motion Diagram for TD1626+368 
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Figure 15. The completeness limit (left) and the proper motion diagram (right) for WD 1626 + 368. 
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Completeness limits for the images of WD1900+705 Proper Motion Diagram for WD1900+705 
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Figure 18. The completeness limit (left) and the proper motion diagram (right) for WD 1900 + 705. 



Searching for Planets Around White Dwarfs 




Completeness limits for the images of WD2007-219 Proper Motion Diagram for WD2007-219 
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Figure 20. The completeness limit (left) and the proper motion diagram (right) for WD 2007 — 219. 
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Completeness limits for the images of WD2140+207 Proper Motion Diagram for WD2140+207 
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Figure 22. The completeness limit (left) and the proper motion diagram (right) for WD 2140 + 207. 



Completeness limits for the images of WD2246 + 223 



Proper Motion Diagram for WD2246 + 223 
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Figure 23. The completeness limit (left) and the proper motion diagram (right) for WD 2246 + 223. 




Figure 24. The completeness limit (left) and the proper motion diagram (right) for WD 2326 + 049. 
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Proper Motion vs Apparent J Magnitude for WD1633+433 
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Figure 25. The motion of the objects between the first epoch 
and second epoch images of WD 1633 + 433. The white dwarf has 
a magnitude of J ~ 14 mag. 



epoch image is required to determine if any of these objects 
are genuine common proper motion companions. 



6.5 WD 2007 - 219 

A single object with a magnitude of J ~ 22 mag appeared 
to have a motion similar to that of WD 2007 — 219 between 
the 2002 first epoch image and 2003 second epoch images 
(Figure 26). This candidate companion was detected with a 
SNR of ~ 9 and ~ 11 in the first epoch and second epoch 
images, respectively. The magnitude of the motion of the 
candidate was ~ 440 mas. In comparison, the rms of the 
magnitude of the motion of all point sources in the field 
with 21.5 < J < 22.5 mag is ~ 130 mas. The motion of 
the candidate was > 3 times larger than this mean error 
and was clearly separated from the other faint objects in 
the field. Therefore, both the magnitude and direction of its 
motion was entirely consistent with that of the white dwarf, 
suggesting that the candidate was a true common proper 
motion companion. The candidate would have a mass of 
7 ± 1 Mjup and a projected physical separation of ~ 980 AU 
if it were confirmed to be a common proper motion compan- 
ion to WD 2007 — 219. A third epoch image was recently ac- 
quired to investigate whether this apparently co-moving can- 
didate was real. Unfortunately, the new observation showed 
that this object is in fact a non-moving background object. 
The NIRI 2008 third epoch image is less noisy than the 
NIRI 2003 second epoch image, which has improved the 
completeness limit and shows the candidate to have some 
structure. Therefore, it is likely to be a galaxy. In summary, 
these observations could have detected a companion with a 
mass of 7 ± 1 Mj up , corresponding to an effective tempera- 
ture of ~ 370 K (Baraffe et al. 2003), between a projected 
physical separation of 55 — 831 AU with a 50% probability 
(Table 3). The motion of this white dwarf between the first 
epoch and third epoch images is large enough to confidently 
state that there are no common proper motion companions 
to WD 2007 - 219 within these limits. 



6.6 WD 2140 + 207 

The mass of this DQ white dwarf used throughout this paper 
(0.49 Mq) was determined by including effects from the car- 
bon present in its atmosphere (Dufour, Bergeron & Fontaine 
2005), while Debes et al. (2006) use a white dwarf mass of 
0.62 Mq (Bergeron, Leggett & Ruiz 2001), which was de- 
rived using a pure helium model. In contrast, a companion 
with a mass of Mj up could have been detected if the 
larger white dwarf mass was used to determine the total age 
of the white dwarf. 

6.7 WD 0148 + 467, WD 1134 + 300, 
WD 1900 + 705, WD 2246 + 223 

The 2003 first epoch images of WD 0148 + 467, WD 1134 + 
300, WD 1900 + 705 and WD 2246 + 223 were degraded by 
severe 60 Hz interference (Table 2), which significantly de- 
creased the completeness limit of these images. It is likely 
that this interference has introduced the large scatter in 
the motions of faint objects between the 2003 first epoch 
and 2004 / 2005 second epoch images (Figures 27, 28, 29 
and 30). This suggests that the error on the motion of 
these faint objects is comparable to the motion of the white 
dwarfs. As a result, multiple objects appear to have mo- 
tions similar to the motion of the white dwarfs. In the case 
of WD 2246 + 223, a single object appears to have a mo- 
tion similar to the motion of WD 2246 + 223. The candi- 
date common proper motion companion is detected with 
a SNR of ~ 5 in both epoch images. If the candidate is 
confirmed to be a common proper motion companion to 
WD 2246 + 223, it would have a mass of 9 ± 1 M Jup and 
a projected physical separation of ~ 840 AU. This corre- 
sponds a projected physical separation of ~ 160 AU around 
the main sequence progenitor, assuming an expansion fac- 
tor of Mms/Mwd ~ 5.1 Af /O.97M ~ 5.3. However, due 
to the presence of 60 Hz interference, the measurement of the 
motion of the candidate between the first epoch and second 
epoch images may be inaccurate. In addition, there are two 
other objects which appear to have the same magnitude of 
motion as the candidate, reducing the probability that the 
candidate is a genuine common proper motion companion. 
Therefore, a third epoch image is required, for all four of 
these white dwarfs, to determine if any of these objects are 
genuine common proper motion companions. 



7 DISCUSSION 

No common proper motion companions within the limits 
given in Table 3 were discovered around 18 of the 23 equa- 
torial and northern hemisphere white dwarfs in the DODO 
survey. Of these, 11 white dwarfs were subject to previous 
searches for substellar companions (Table 4). The DODO 
survey extends out to a much larger projected physical sep- 
aration in each case and places a new lower, upper limit on 
the mass of any possible companion in orbit around 7 of 
these white dwarfs (Tables 3 and 4). For the remaining 5 
targets, multiple objects in each field of view appear to have 
motions similar to the motion of the white dwarfs. 
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Table 3. Results for the 23 equatorial and northern hemisphere white dwarfs in the DODO survey 



White 


hot 


90% J 


90% M 


90% T 


50% J 


50% M 


50% T 


WD Orbit 


MS Orbit 


lcr error 




Dwarf 


[Gyr] 


[mag] 


[M Jup ] 


[K] 


[mag] 


[M Jup ] 


[K] 


[AU] 


[AU] 


[mas] 




0115 


+ 159 


1.7 


21.0 


10+? 
lu -i 


430 


22.0 


8+} 
°— l 


380 


46 - 


675 


11 - 


160 


94 




0148 


+ 467 


2.5 


20.4 


16+| 


180 


21.9 


10+ 2 


390 


48 - 


457 


14 - 


138 


66 


1.2 


0208 


+ 396 


2.6 


20.5 


16+ 3 


480 


22.5 


9+} 


360 


50 - 


758 


13 - 


194 


91 




0341 


+ 182 


3.3 


22.0 


13+3 


400 


22.9 


10+ 2 

- 1 


360 


57 - 


801 


16 - 


222 


79 




0435 


-088 


4.1 


21.2 


13+3 


380 


22.7 


9+i 


320 


28 - 


408 


9 - 


124 


94 




0644 


+ 375 


2.1 


20.5 


13+2 


460 


22.4 


8 +1 
°-i 


360 


46 - 


652 


17 - 


236 


124 


1 


0738 


- 172 


2.4 








22.0 


7+1 


320 


27 - 


379 


7- 


96 


42 


3 


0912 


+ 536 


3.0 


20.9 


13+° 


110 


22.1 


9+i 


350 


31 - 


419 


7 - 


93 


171 


1 


1055 


-072 


3.3 


21.0 


13+3 


400 


22.6 


9 +1 
— l 


340 


36 - 


503 


8 - 


103 


80 




1121 


+ 216 


2.3 


21.2 


10+ 2 
lu -i 


390 


22.2 


g+2 
°- 1 


350 


40 - 


605 


9 - 


134 


117 


1 


1134 


+ 300 


0.37 


20.8 


5 +0 
J -i 


440 


21.9 


3+i 

°-0 


350 


46 - 


664 


9 - 


127 


97 


1.2 


1344 


+ 106 


2.5 


20.8 


16+1 


480 


22.0 


13+o 


440 


60 - 


865 


14 - 


208 


111 


1 


1609 


+ 135 


2.8 


21.7 


13±i 


420 


22.5 


10+ 2 


380 


55 - 


642 


10 - 


117 


83 


1 


1626 


+ 368 


2.2 


22.1 


9±\ 


380 


22.8 




360 


48 - 


535 


13 - 


141 


72 


1 


1633 


+ 433 


3.0 


21.1 


13±g 


410 


22.3 


10+1 


370 


45 - 


533 


10 - 


123 


97 


1,2,4 


1647 


+ 591 


0.91 


19.6 


9±i 


480 


22.0 


5+? 


350 


33 - 


372 


7- 


77 


127 




1900 


+ 705 


1.1 


21.2 


Tt\ 


100 


22.2 


5±S 


330 


39 - 


452 


8 - 


89 


146 


1.2 


1953 


-011 


2.1 


19.2 


1615 


510 


21.7 




360 


34 - 


509 


7 - 


111 


60 


1 


2007 


-219 


1.4 


21.2 




450 


22.4 


rt\ 


370 


55 - 


831 


12 - 


189 


74 




2047 


+ 372 


0.89 








21.8 


6+ 1 
°-0 


390 


54 - 


202 


12 - 


46 


11 


1 


2140 


+ 207 


4.4 


20.0 


21+° 


490 


21.6 


13±g 


370 


38 - 


542 


13 - 


181 


73 


1 


2246 


+ 223 


1.7 


20.6 


13±g 


490 


22.0 


»±1 


400 


57 - 


835 


11 - 


157 


92 


1.2 


2326 


+ 049 


1.1 


21.1 


7±\ 


400 


21.8 


«tl 


370 


41 - 


396 


9 - 


89 


110 


1.5 



Columns: ttot is the "COND" evolutionary model age used; 90% and 50% gives the 90% and 50% completeness limits in terms 
of apparent J magnitude, mass, M, measured in Jupiter masses, and effective temperature, T, measured in Kelvin, respectively; 
WD Orbit is the range of projected physical separations at which a companion of that mass could be found around the white 
dwarf, measured in AU; MS Orbit is the range of projected physical separations at which a companion of that mass could be found 
around the main sequence progenitor, measured in AU; lcr error is the lcr scatter of the distribution of the motions of all objects 
in the field, excluding the white dwarf, measured in milli arc seconds. Notes: (1) 60Hz signal present in at least one epoch image, 
(2) A third epoch image is required, (3) A common proper motion companion is known to orbit this white dwarf, (4) Streak across 
one epoch image due to a bright source just outside the field of view, (5) A dust disk is known to orbit this white dwarf. 



Proper Motion Diagram for WD2007-Z19 



Proper Motion Diagram for WD2007-219 
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Figure 26. The proper motion diagram (left) shows the motion of all objects in the field of WD 2007 — 219 between the first epoch and 
second epoch images. The dashed green circles represent the lcr and 3a scatter of the distribution of the motions of all objects excluding 
the white dwarf, centred on the white dwarf, to help determine possible common proper motion companions to the white dwarf. A single 
object appears to have a motion similar to that of WD 2007 — 219. The latest proper motion diagram (right), utilising the 2008 third 
epoch image, shows that there are no common proper motion companions to WD 2007 — 219. 
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Table 4. Previous searches for substcllar companions around the 
23 equatorial and northern hemisphere white dwarfs in the DODO 
survey 



Proper Motion vs Apparent J Magnitude for WDO 148+467 



White 


M 


WD Orbit 


WD Orbit 


Ref 


Dwarf 


[M Jup ] 


["] 


[AU] 




01 1 ^ -U 1 










uiio '-tu i 










0908 -1- 3Qfi 


10 


0.9 - 10 


15 - 167 


i 


0^41 _L 1 89 










OzL^ — 088 
yioo uoo 










0fi4zL -U 37^ 










07^8 — 1 79 
U 1 oo J. / ^ 










0Q1 9 -U ^fi 


12 


1 - 7 


10 - 72 


2 


10^ — 079 


14 


1 - 7 


12 - 85 


2 


1121 + 216 


11 


1 - 7 


13 - 94 


2 


1134 + 300 










1344 + 106 


14 


Unresolved 


3 


1609 + 135 










1626 + 368 


14 


1 - 7 


16 - 112 


2 


1633 + 433 


14 


1 - 7 


15 - 106 


2 




14 


Unresolved 


3 


1647 + 591 










1900 + 705 










1953-011 


10 


1 - 7 


11 - 80 


2 


2007 - 219 










2047 + 372 










2140 + 207 


10 


1 - 7 


13 - 88 


2 


2246 + 223 


9 


1 - 7 


19 - 133 


2 


2326 + 049 


6 


1 - 5 


14 - 68 


4 



Columns: M is the minimum mass of a companion that could 
be found around the white dwarf, measured in Jupiter masses; 
WD Orbit is the range of projected physical separations at 
which a companion of that mass could be found around the 
white dwarf, measured in arc seconds and determined in AU 
using the distance to the white dwarf, measured in parsccs, 
taken from van Altcna, Lee & HofHeit (1995) (Table 1); Two 
targets were part of a search for companions through the de- 
tection of an infrared excess and are denoted by "Unresolved" ; 
Ref = References: (1) Debes, Sigurdsson & Woodgatc (2005b), 
(2) Debes, Ge & Ftaclas (2006), (3) Farihi, Bccklin & Zuck- 
crman (2008), (4) Debes, Sigurdsson & Woodgate (2005a). 
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Apparent J Magnitude (mag) 

Figure 27. The motion of the objects between the 2003 first 
epoch and 2005 second epoch images of WD 0148 + 467. The 
white dwarf has a magnitude of J ~ 13 mag. 



Proper Motion vs Apparent J Magnitude for WD1134+300 



Apparent J Magnitude (mag) 

Figure 28. The motion of the objects between the first epoch 
and second epoch images of WD 1134 + 300. The white dwarf has 
a magnitude of J ~ 13 mag. 



Two of these stars, WD 0148 + 467 and WD 1134 + 
300, have the lowest proper motions in the survey, moving 
only ~ 225 mas and ~ 286 mas, respectively, between their 
2003 first epoch and 2005 second epoch images. This motion 
is only ~ 1.7 and ~ 2.0 times larger than the rms of the 
magnitude of the motion of all the point sources in the fields 
of WD 0148 + 467 and WD 1134 + 300, respectively, with 
21.5 < J < 22.5 mag. In addition, their first epoch images 
were degraded by 60 Hz interference. Due to the combination 
of the very low proper motions of these white dwarfs and 
the presence of 60 Hz interference, third epoch images are 
required to clearly distinguish their proper motions from the 
motions of the background objects in the field. 

WD 1633 + 433, WD 1900 + 705 and WD 2246 + 223 
all move > 450 mas between their 2003 first epoch and 2004 
second epoch images. However, their first epoch images were 
degraded by 60 Hz interference. In addition, the presence of 
a streak in the first epoch image of WD 1633 + 433 has 
most likely decreased the accuracy of the measurement of 
the motion of the faintest objects in the field. Again, third 



epoch images, without the presence of any interference, are 
required. 

Using the 2008 third epoch image of WD 2007 - 219, 
the necessary minimum baseline to reliably distinguish real 
common proper motion companions from non-moving back- 
ground objects was determined. For the faintest (and by 
implication, the most interesting) sources (SNR < 10, cor- 
responding to J > 22 for the majority of the images in 
the DODO survey) the white dwarf needs to have moved 
at least 4 pixels (~ 470 mas for the NIRI data) between 
the two epoch images to conclusively rule out non-moving 
outliers at this magnitude. 

The cumulative completeness limits, in terms of mass 
and effective temperature, and the corresponding range of 
projected physical separations over which these limits ap- 
ply have been determined for all 23 equatorial and north- 
ern hemisphere white dwarfs discussed in this paper (Fig- 
ures 31, 32, 33, 34). 

From these results, tentative conclusions regarding the 
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Proper Motion vs Apparent J Magnitude for WD1900+705 




Apparent J Magnitude (mag) 

Figure 29. The motion of the objects between the first epoch 
and second epoch images of WD f 900 + 705. The white dwarf has 
a magnitude of J ~ 13.5 mag. 



Proper Motion vs Apparent J Magnitude for WD2246+223 



Apparent J Magnitude (mag) 



Figure 30. The motion of the objects between the first epoch 
and second epoch images of WD 2246 + 223. The white dwarf has 
a magnitude of J ~ 14.3 mag. 



frequency of substellar and planetary mass companions to 
white dwarfs and their main sequence progenitors at wide 
separations can be made (we recognise that the DODO sur- 
vey contains a relatively small number of targets). These 
conclusions assume that no common proper motion com- 
panions are confirmed around the 5 white dwarfs requiring 
a third epoch image and include the non-detection of a com- 
panion around WD 0046 + 051 (Burleigh et al. 2008). Firstly, 
using the 90% completeness limits, the DODO survey can 
detect companions with effective temperatures > 500 K 
around all targets. This is significantly below the currently 
coolest known brown dwarfs, ULAS J003402.77 - 005206.7 
(Warren et al. 2007) and CFBDS J005910.90 - 011401.3 
(Delorme et al. 2008), which have effective temperatures 
of 600 < T eff < 700 K and spectral types of T8.5. In 
fact, these observations probe well into the hypothetical 
Y dwarf regime or at least into a significant extension of 
the T dwarf sequence (Kirkpatrick 2005). Therefore, we 
suggest that < 5% of white dwarfs have L, T and sub- 



T8.5 (Y?) substellar companions with effective tempera- 
tures > 500 K between projected physical separations of 
60 — 200 AU, although for many fields this applies to smaller 
(~ 13 AU for WD 0046+051; Burleigh et al. 2008) and larger 
(~ 800 AU) projected physical separations. This corre- 
sponds to projected physical separations around their main 
sequence progenitors (1.5 — 8 M©, i.e., spectral types F5-B5) 
of 20 — 45 AU, although again for many fields these limits 
apply to smaller (~ 3 AU for WD 0046 + 051; Burleigh et al. 
2008) and larger (~ 200 AU) projected physical separations. 
For the same range of projected physical separations stated 
above and using the 50% completeness limits, we suggest 
that < 8% of white dwarfs and their main sequence pro- 
genitors have companions with masses above the deuterium 
burning limit (~ 13Mj up ), while < 9% have companions 
with masses > 10Mj up . 

These results can be compared to the results from other 
imaging surveys for wide substellar and planetary mass com- 
panions to white dwarfs and main sequence stars (Table 5). 
In particular, our results are consistent with those of Mc- 
Carthy & Zuckerman (2004) and Lafreniere et al. (2007). We 
note the recent claims of the directly imaged planetary mass 
companions to Fomalhaut (Kalas et al. 2008) and HR8799 
(Marois et al. 2008) and await the statistical analyses of 
those surveys for comparison with the DODO survey re- 
sults. The DODO survey results can also be compared to 
complimentary recent MIR searches for unresolved substel- 
lar and planetary mass companions to white dwarfs (e.g., 
Mullally et al. 2007). A recent MIR photometric survey of 
27 white dwarfs using the Spitzer Space Telescope and IRAC 
was sensitive to the entire known T dwarf sequence (Farihi 
et al. 2008). Their observations place similar limits (< 4%) 
on the frequency of such companions to white dwarfs, but 
at smaller separations (with some overlap) compared to the 
DODO survey. 

At this stage, we prefer to refrain from speculating on 
the reasons for the negative results so far. For example, from 
radial velocity measurements of evolved giant stars, Lovis & 
Mayor (2007) estimate that at least 3% of stars with M > 
1.8 Mq host M p sin i > 5 M.j up companions, including brown 
dwarfs. Therefore, it is not unreasonable to have expected 
that at least one of our target white dwarfs would have had a 
detectable companion. It may simply be a case of observing 
more targets. 

A more substantial comparison with similar surveys, 
along with a more thorough statistical analysis of our results, 
will be presented in a forthcoming paper on the southern 
hemisphere white dwarfs in the DODO survey. 
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Table 5. Recent imaging searches for wide companions 



Survey 


Targets 


Number 
of Targets 


Limit 
[M Jup ] 


Separation 
[AU] 


Frequency of 
Companions 


McCarthy & Zuckerman (2004) 


G K M 


102 


> 12 


75 - 


300 


1% ± 1% 






178 


> 30 


140 - 


1200 


0.7% ±0.7% 








5-10 


75 - 


300 


< 3% 


Farihi et al. (2005) 


White Dwarfs 


261 


> 52 


100 - 


5000 


< 0.5% 






86 


> 21 


50 - 


1100 


< 0.5% 


Allen et al. (2007) 


M7-L8 


132 


> 52 


40 - 


1000 


< 2.3% 


Lafreniere et al. (2007) 


F G K M 


85 


13-40 


25 - 


250 


< 5.6% 


Nielsen et al. (2008) 


A F G K M 


(50 


> 4 


20 - 


100 


< 20% 
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20 
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Figure 31. The cumulative completeness limit, in terms of com- 
panion mass, for the 23 equatorial and northern hemisphere white 
dwarfs in the DODO survey. The red dotted-dashed line indicates 
the frequency of the completeness limit in M, j up at which 90% of 
companions with that mass could be detected, while the blue 
dashed line indicates the completeness limit in Mj up at which 
50% of companions with that mass could be detected. 
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Figure 32. The cumulative completeness limit, in terms of com- 
panion temperature, for the 23 equatorial and northern hemi- 
sphere white dwarfs in the DODO survey. The red dotted-dashed 
line indicates the frequency of the completeness limit in Kelvin 
at which 90% of companions with that temperature could be de- 
tected, while the blue dashed line indicates the completeness limit 
in Kelvin at which 50% of companions with that temperature 
could be detected. 
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Figure 33. The cumulative minimum projected physical separa- 
tions for the 23 equatorial and northern hemisphere white dwarfs 
in the DODO survey. The red dotted— dashed line indicates the 
minimum projected physical separations in AU at which a com- 
panion could be found around each white dwarf. The blue dashed 
line indicates the minimum projected physical separations in AU 
at which a companion could be found around the main sequence 
progenitor. 
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Figure 34. The cumulative maximum projected physical separa- 
tions for the 23 equatorial and northern hemisphere white dwarfs 
in the DODO survey. The red dotted— dashed line indicates the 
maximum projected physical separations in AU at which a com- 
panion could be found around each white dwarf. The blue dashed 
line indicates the maximum projected physical separations in AU 
at which a companion could be found around the main sequence 
progenitor. 
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